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ABSTRACT 

We present a study of the spatial distribution of dwarf satellites (or subhalos) in galactic dark matter halos using 
dissipationless cosmological simulations of the concordance flat Cold Dark Matter (CDM) model with vacuum 
energy. We find that subhalos are distributed anisotropically and are preferentially located along the major axes 
of the triaxial mass distributions of their hosts. The Kolmogorov-Smirnov probability for drawing our simulated 
subhalo sample from an isotropic distribution is Pks — 1.5 x 10~ 4 . An isotropic distribution of subhalos is thus 
not the correct null hypothesis for testing the CDM paradigm. The nearly planar distribution of observed Milky 
Way (MW) satellites is marginally consistent (probability ~ 0.02) with being drawn randomly from the subhalo 
distribution in our simulations. Furthermore, if we select the subhalos likely to be luminous, we find a distribution 
that is consistent with the observed MW satellites. In fact, we show that subsamples of the subhalo population with 
a centrally-concentrated radial distribution that is similar to that of the MW dwarfs typically exhibit a comparable 
degree of planarity. We explore the origin of the observed subhalo anisotropy and conclude that it is likely due 
to (1) the preferential accretion of satellites along filaments, often closely aligned with the major axis of the host 
halo, and (2) evolution of satellite orbits within the prolate, triaxial potentials typical of CDM halos. Agreement 
between predictions and observations requires the major axis of the outer dark matter halo of the Milky Way to be 
nearly perpendicular to the disk. We discuss possible observational tests of such disk-halo alignment with current 
large galaxy surveys. 

Subject headings: cosmology: theory, large-scale structure of universe - dark matter - galaxies: formation, 
halos, structure — methods: numerical 



1. INTRODUCTION 

Simulations of structure formation in the standard cold dark 
matter (CDM) scenario (e.g., Blumenthal et al. 1984) show that 
virialized dark matter halos teem with distinct, gravitationally- 
bound substructures, often referred to as subhalos. The abun- 
dance of subhalos in Milky Way-sized halos has received much 
attention, as there are more than an order of magnitude fewer 
observed dwarf satellite galaxies around these systems than the 
predicted number of subhalos of comparable velocity disper- 
sion (Kauffmann et al. 1993; Klypin et al. 1999a; Moore et al. 
1999). This problem is often referred to as the "the missing 
satellites problem." 

The difference in predicted and observed abundances of 
dwarf satellites is likely related to the corresponding differ- 
ence in spatial distributions. The predicted radial distribu- 
tion of subhalo populations is considerably more extended than 
that of the Milky Way (MW) satellites (Taylor et al. 2003; 
Kravtsov et al. 2004). This may indicate that the environments 
in which low-luminosity galaxies form are non-trivially biased 
relative to the overall population of subhalos. In addition, the 
dwarf satellites of the MW and M31 appear to be distributed 
anisotropically about their hosts (e.g., Lynden-Bell 1982; Ma- 
jewski 1994; Hartwick 1996, 2000; Mateo 1998; Grebel et al. 
1999; Willman et al. 2004). Kroupa et al. (2005) recently ar- 
gued that the anisotropic distribution of the MW dwarf satel- 
lites presents a serious challenge to the standard CDM structure 
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formation paradigm. Kroupa et al. (2005) reached this conclu- 
sion by assuming that the luminous satellites correspond to a 
randomly-selected subset of dark matter subhalos and by tak- 
ing an isotropic subhalo distribution as their null hypothesis. 

There have been a number of studies of anisotropy in satel- 
lite galaxy distributions outside of the Local Group. Holm- 
berg (1969) found that satellites of spiral galaxies with pro- 
jected separations r p < 50kpc are preferentially located near 
the short axes of the projected light distributions of their host 
galaxies (the Holmberg Effect). However, several subsequent 
studies found little evidence for such a preferential alignment 
(e.g., Hawley & Peebles 1975; Sharp et al. 1979; MacGillivray 
et al. 1982). Zaritsky et al. (1997) found a statistically signifi- 
cant anisotropy similar to that advocated by Holmberg (1969), 
but only at larger projected distances, 200kpc< r p < 500 kpc. 
In a more recent study of a large sample of satellite galax- 
ies in the Two Degree Field Galaxy Redshift Survey (Colless 
et al. 2001), Sales & Lambas (2004) also found evidence for 
preferential alignment of satellites along the minor axis of the 
central galaxy. However, Brainerd (2004) performed a simi- 
lar study of satellites in the Sloan Digital Sky Survey (York 
et al. 2000; Strauss et al. 2002), finding satellites to be pref- 
erentially aligned near the long axes of host galaxies and the 
degree of anisotropy to be a rapidly decreasing function of sep- 
aration from the host galaxy. Note that the anisotropy measured 
by Brainerd (2004) is opposite in sense to the anisotropy re- 
ported by Holmberg (1969), Zaritsky et al. (1997), and Sales & 
Lambas (2004). 

Properties of dwarf satellite dark matter (DM) halos (or sub- 
halos) in MW-sized hosts have been the subject of several re- 
cent studies, which used a new generation of high-resolution 
dissipationless simulations not affected by the "overmerging" 
problem (e.g., Klypin et al. 1999a; Moore et al. 1999; Stoehr 
et al. 2002, 2003; De Lucia et al. 2004; Kravtsov et al. 2004; 
Gao et al. 2004b; Reed et al. 2004). One of the main results 
is that the radial distribution of subhalos is more extended than 
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the distribution of DM (Ghigna et al. 1998; Colin et al. 1999; 
Ghigna et al. 2000; Springel et al. 2001; De Lucia et al. 2004; 
Gao et al. 2004a). The DM subhalos also appear to have a sig- 
nificantly more extended and shallower radial distribution com- 
pared to the observed distribution of satellite galaxies both in 
galactic halos (Taylor et al. 2003; Kravtsov et al. 2004) and in 
cluster halos (Diemand et al. 2004; Gao et al. 2004a; Nagai & 
Kravtsov 2005). Despite significant progress in our understand- 
ing of halo substructure, the anisotropy of subhalo spatial dis- 
tributions has, so far, not been studied in as much detail. Zarit- 
sky et al. (1997) reported no detectable anisotropy in the pro- 
jected satellite distribution in the simulations of Navarro, Frenk, 
& White (1994, 1995). However, this conclusion refers to the 
statistical, projected distribution of the most massive satellites 
obtained after stacking many galaxy-sized halos, rather than 
anisotropy of the satellite distribution within a single host halo. 
More recently, Knebe et al. (2004) studied anisotropy of the 
subhalo distribution in dissipationless simulations of cluster- 
sized hosts. These authors found that subhalo distribution is 
anisotropic and is aligned with the major axis of the matter dis- 
tribution of the host. 

In this paper, we study the anisotropy of satellite distribu- 
tion in Galaxy-sized halos using high-resolution cosmological 
iV-body simulations of structure formation in the concordance 
flat ACDM cosmology. We show that an isotropic distribu- 
tion is not the correct null hypothesis for testing the CDM 
paradigm. The mass distributions in CDM halos are gener- 
ally triaxial rather than spherical. We demonstrate that sub- 
halos of the size thought to host the MW dwarf satellites are 
distributed anisotropically about their host halos with subhalos 
preferentially located along the major axes of their hosts. We 
also show that the null hypothesis distribution taken by Kroupa 
et al. (2005) is incorrect even in the case of an isotropic un- 
derlying distribution. As we were completing this study, Kang 
et al. (2005) presented a similar study, considering in particu- 
lar the question of whether the anisotropy of the MW satellite 
distribution is consistent with the hierarchical formation sce- 
nario. Although our approaches differ in detail, their conclu- 
sions are consistent with ours. We present a more extensive 
study of the satellite distribution, both in three dimensions and 
in two-dimensional projection. We also explore the physical 
mechanisms that create the anisotropy measured in the simula- 
tions and discuss the implications of our results on our under- 
standing of galaxy formation. 

This manuscript is organized as follows. In § 2, we describe 
our numerical simulations and our analysis methods. In § 3, 
we present results on the anisotropic distribution of subhalos in 
Galaxy-sized halos. We discuss our results and their implica- 
tions in § 4. Lastly, we summarize our main findings and draw 
conclusions in § 5. In an Appendix, we discuss the prospects of 
studying satellite anisotropy in projection. 

Throughout this paper, we refer to halos that are contained 
within the virial radii of still larger halos as subhalos or satel- 
lites and we refer to halos that are not contained within a larger 
halo as host halos. 

2. METHODS 

2.1. Numerical Simulations 

We analyze a simulation of three Milky Way-sized DM ha- 
los formed in a standard, "concordance" ACDM cosmology 
with n M =l-n A = 0.3, h = 0.7, and cr 8 = 0.9. The simula- 
tion was performed with the Adaptive Refinement Tree (ART) 



Af-body code (Kravtsov et al. 1997; Kravtsov 1999). This sim- 
ulation has been discussed previously by Klypin et al. (2001), 
and Kravtsov, Gnedin, & Klypin (2004, hereafter KGK04). We 
briefly review some of the details in this section. 

The simulation began with a uniform 256 3 grid covering a 
comoving, cubic box of 25 /z _1 Mpc on a side. Higher force 
resolution was achieved in dense regions using an adaptive re- 
finement algorithm. The grid cells were refined if the parti- 
cle mass contained within them exceeded a certain, specified 
threshold value. The multiple mass resolution technique was 
used to set up the initial conditions. A Lagrangian region corre- 
sponding to two virial radii about each halo was re-sampled at 
the initial epoch of z, = 50 with the highest-resolution particles 
ofmassm p = 1.2 x lO 6 /z _1 M . The high mass resolution region 
was surrounded by layers of particles of increasing mass with 
a total of five particle species. Only the regions containing the 
highest-resolution particles were adaptively refined. The max- 
imum level of refinement in the simulations corresponded to a 
peak formal spatial resolution of approximately 100/i _1 pc. 

We define a halo radius /?iso, by the sphere, centered on the 
particle with the highest density, within which the mean density 
is 180 times the mean density of the universe p M , so that the 
mass and radius are related by M\%q = 47t(180/9m)^i8o/3- Our 
three host halos, which we refer to as halo G\, halo G2, and 
halo G3, have masses of M m) = 1.66 x 10 12 /t'Mq, 1.24 x 
10 12 h~ l M Q , and 1.19 x 10 12 h~ l M Q , respectively, and these 
halos contain ~ 10 6 particles within their virial radii. 5 Their 
virial radii are R m = 298/r'kpc, 278/T 1 kpc, and 281/r'kpc, 
respectively. 

We identified halos and subhalos using a variant of the Bound 
Density Maxima algorithm (Klypin et al. 1999b). First, we 
compute the local density at each particle using a smoothing 
kernel of 24 particles and identify local maxima in the density 
field. Beginning with the highest density particles and stepping 
down in density, we mark each peak as a potential halo center 
and surround the peak by a sphere of radius r^ = 10/z _1 kpc. 
All particles within the sphere are excluded from further con- 
sideration as potential halo centers. The parameter r^ is set 
according to the smallest objects that we aim to identify ro- 
bustly. After identifying potential halo centers, we iteratively 
remove unbound particles. For host halos, the mass and radius 
are set according to a fixed overdensity as described above. For 
subhalos, the outer boundary is somewhat ambiguous and we 
adopt a truncation radius r t , at which the density profile be- 
comes greater than a critical value of dlnp/dlnr = -0.5. This 
criterion is based on the fact that we do not expect density pro- 
files of CDM halos to be shallower than this and, empirically, 
this definition is approximately equal to the radius at which the 
background density of host halo particles is equal to the density 
of particles bound to the subhalo. 

Upon identifying halos, we assign each halo a mass and ra- 
dius and use the halo particles to determine a circular velocity 
profile V c i rc (r) = y/GM{< r)/r, and the maximum circular ve- 
locity, V max - We choose to quantify the size of subhalos accord- 
ing to y max because this quantity is measured more robustly and 
is not subject to the same ambiguity as a particular mass defini- 
tion. 

Halos Gi, G2, and G3 have maximum circular velocities of 
V max ~ 213kms- 1 , V max ~ 199 kms" 1 , and V max ~ ^kms" 1 re- 
spectively. These halos were selected to reside in a well-defined 

5 These halos were referred to as B2, C2, and D2 respectively by Klypin et al. 
(2001). 
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filament at z = 0. The halos G2 and G3 are neighbors located at 
425/T 1 kpc (i.e., w 610 kpc ~ 2/? v i r ) from each other. The con- 
figuration of this pair thus resembles that of the Local Group. 
The third halo is isolated and is located ~ 2 Mpc away from 
the pair. The three hosts have similar masses at the present but 
rather different mass accretion histories. Host Gi undergoes a 
spectacular multiple major merger at z ~ 2, which results in a 
dramatic mass increase on a dynamical time scale. Halos G2 
and G3 increase their mass in a series of somewhat less spec- 
tacular major mergers which could be seen as mass jumps at 

5 < z < 1 . All three systems accrete little mass and experience 
no major mergers at z < 1 (or lookback time of w 8 Gyr) and 
thus could host a disk galaxy. Note, however, that halos Gi and 
G3 experience minor mergers during this period. 

2.2. Modeling Luminous Satellites 

The relative scarcity of MW satellites compared to predicted 
subhalo counts suggests that a naive comparison of subhalo 
populations to luminous galaxies may not be correct. This im- 
plies that in order to make a more meaningful comparison with 
the data, we should have a model to identify the DM subhalos 
in simulations that would host observable, luminous galaxies. 
We consider two alternatives. 

The first dwarf galaxy formation model we consider was re- 
cently proposed by KGK04 (see their § 6 for details). This 
semi-analytic model is based on the subhalo evolutionary tracks 
extracted from the simulations used in this study. The small- 
mass dwarfs are identified with the halos that either have ac- 
creted a large fraction of their mass prior to the epoch of reion- 
ization (see also Bullock, Kravtsov, & Weinberg 2000; Ricotti 

6 Gnedin 2004) or were relatively massive at high redshifts 
and could therefore retain most of their gas and form stars af- 
ter reionization. Some objects in the latter category could lose 
most of their former mass due to tidal stripping and appear as 
relatively low-mass halos at the present epoch. The dwarf satel- 
lite galaxies in this model can thus be hosted by both massive 
and low-mass subhalos at z = 0. 

The physical ingredients of the KGK04 model include: (1) 
the suppression of gas accretion by the extragalactic UV back- 
ground; (2) an observationally-motivated recipe for quiescent 
star formation; (3) a starburst mode of star formation af- 
ter strong peaks in tidal forces (which are calculated self- 
consistently from the simulations); and (4) an accounting for 
the inefficient dissipation of gas in halos with T^ n < 10 4 K. The 
model can successfully reproduce the circular velocity function, 
the radial distribution, and the morphological segregation of the 
observed MW satellites, as well as the basic properties of galac- 
tic dwarfs such as their star formation histories, stellar masses, 
and densities. In our analysis, we use the same set of subhalos 
that were considered to be "luminous" according to this model 
in KGK04 and refer to them as luminous subhalos. 

The second model assumes that the observed MW satellites 
are hosted by the most massive subhalos (Stoehr et al. 2002, 
2003). In this model, it is conjectured that the masses of the 
subhalos in which the luminous dwarfs are embedded are sig- 
nificantly underestimated because the DM density profiles in 
the central regions of the subhalos have been affected by tidal 
interactions (Stoehr et al. 2002, 2003; Hayashi et al. 2003, see, 
however, Kazantzidis et al. 2004b). The maximum circular ve- 
locities of the Local Group dwarfs may thus be systematically 
underestimated in observations because they are derived from 
stellar velocity dispersion measurements within radii consider- 
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FIG. 1 . — Host halo axis ratio profiles for the three MW-sized host halos. 
Top: intermediate-to-long axis ratio q = b/a, as a function of long axis length. 
Bottom: short-to-long axis ratio s = c/a. In both panels, the solid line rep- 
resents halo Gi, the dashed line represents halo G2, and the dot-dashed line 
represents halo G3 . The thick lines represent the shape profiles at z = 0, while 
the thin lines represent the halo shape profiles at z = 1 . 



ably smaller than the radius at which V max is achieved. Stoehr 
et al. (2002) argued that the bias induced by this is large, such 
that all of the observed MW satellites can be embedded in the 
most massive subhalos with V mdx > 30-40kms _1 . This model 
has an important physical implication. If the MW dwarfs in- 
habit the most massive subhalos, then there must exist a certain 
universal mass or V max scale below which galaxy formation is 
completely quashed due to the UV background heating and in- 
efficient gas cooling in dwarf halos. We consider this type of 
model by associating luminous dwarfs with the eleven subhalos 
with the highest values of V max at z = within 300 kpc of each 
MW-sized host halo. 

2.3. The Principal Axes of the Host Halos 

We determine the principal axes of the three simulated host 
halos and the corresponding principal axis ratios q = b/a and 
s = c/a (a > b > c) in the following way. We construct a mod- 
ified inertia tensor given by (e.g., Dubinski & Carlberg 1991) 

7y=$>„*r*y/& (i) 

V 

where m v is the mass of the i4h particle, x\ is the i coordinate 
with respect to a reference frame centered on the density peak 
of the halo, Cl = (y\) 2 + (y% /s) 2 + (y% /q) 2 , and y\ are the parti- 
cle coordinates with respect to the halo principal axes. We use 
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FIG. 2. — Projections of satellites in a plane orthogonal to their best-fit planes (see text). All panels show scatter plots of the positions of satellites projected onto 
a plane perpendicular to their best-fit plane. The best-fit plane corresponds to the vertical axis in this projection. In both rows, the first panel shows results for the 
subhalos of halo G\, the second column for halo G2, and the third column for halo G3. The fourth column in each row shows the observed MW satellites. In the 
first three columns, the projection is such that the major axis of the host halo lies in the plane of the projection. In these panels, the major axis is shown as the thin, 
solid line. In the fourth column, the projection is such that the MW disk is seen edge-on and the MW disk orientation is denoted by the thick, solid line. In the top 
row, we compute the best-fit plane by considering all subhalos within sa 300 kpc of the center of the host halo. In the bottom row, we compute the best-fit plane with 
respect to all luminous subhalos within 30Qkpc of the host halo center. 



an iterative algorithm to determine the principal axes. We be- 
gin with the assumption of a spherical configuration (a = b = c), 
construct the inertia tensor according to Equation (1), and di- 
agonalize the tensor to determine the principal axes (eigenvec- 
tors) and the axis ratios (ratios of eigenvalues). We then repeat 
this process, using the results of the previous iteration to de- 
fine the principal axes, until the results converge to a fractional 
difference of 10~ 3 . The process generally takes fewer than 10 
iterations to converge. The factor of £~ 2 in Eq. (1) serves to mit- 
igate the influence of massive substructures at large distances, 
which can be a significant source of noise in the measurement 
(Dubinski & Carlberg 1991). 

In Figure 1, we show shape profiles at z = and z = 1 for each 
host halo constructed in this way. The Figure shows the axis ra- 
tios as a function of the length of the major axis of the ellipsoid 
£, and we construct profiles by considering all particles within 
this ellipsoid. In agreement with previous studies (e.g., Jing & 
Suto 2002; Bullock 2002), all three MW-sized host halos are 
triaxial with b/a k 0.6-0.7 and tend to be more prolate than 
oblate with c/b rs 0.8-0.95. Additionally, it is evident that halo 
shapes evolve as halos tend to be less spherical at high redshift, 
reflecting their younger dynamical age. 

Currently, it is unclear whether such axis ratios are consistent 
with observational constraints on the shape of the MW halo. 
Recent studies suggest that the coherence of the Sagittarius tidal 
debris constrains the minor- to-major ratio to c/a> 0.8 (Ibata 
et al. 2001; Majewski et al. 2003; Johnston et al. 2004). How- 
ever, Helmi (2004a,b) argues that the Sagittarius tidal stream is 
consistent with a minor- to-major axis ratio as small as c/a ~ 
0.6 (see also Martfnez-Delgado et al. 2004) because the stream 
is dynamically young. Moreover, the shape constraint could be 



relaxed in a prolate halo where orbits along the long-axis are 
less susceptible to strong precession (see, however, Johnston 
et al. 2004). We note that the effect of baryon cooling during 
galaxy formation should make halos more spherical compared 
to halos in the dissipationless simulations analyzed here. The 
change in axis ratios is A(c/a) ~ 0.1-0.3 in the inner regions 
of the halo, and it is not uncommon to see significant variation 
in axis ratios with ( (Dubinski 1994; Kazantzidis et al. 2004a). 

In our analysis below, we refer to subhalo positions in a coor- 
dinate system defined by the principal axes of the host halo in- 
ertia tensor, calculated using DM particles within the ellipsoid 
with ( = 0.3/?i8o (for the MW-sized halos, this corresponds to 
£ w 85/r'kpc w 120 kpc). We choose particles with £ < 0.3/?iso 
in order to mitigate the influence of large substructures at large 
halo-centric distances on the reference frame definition; how- 
ever, we find that the ellipsoids defined at different values of £ 
are closely aligned at all radii within any single host, in agree- 
ment with Jing & Suto (2002). We define the zenith angle, 
< 9 < 7r/2, as the angle from the major axis of the halo. We 
also make analogous computations in two-dimensional projec- 
tion, in which case we define a polar coordinate system with 
angle <j), defined as the angle away from the major axis of the 
two-dimensional projected DM distribution. 

2.4. Satellite Planes 

Instead of defining satellite positions in the coordinate sys- 
tem set by the mass distribution of their host halo, which is 
difficult to determine observationally, one can construct coordi- 
nate systems with respect to the satellites themselves. Kroupa 
et al. (2005) found the MW satellites to be in a nearly planar 
distribution, so we follow the method of Kroupa et al. (2005) 
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FIG. 3. — The cumulative fraction of halos with an angular position < | cos(cj)| with respect to the normal to the best-fit plane as a function of | cos(w)|. The left 
panel shows the distribution for all subhalos with V max > 12kms~', the center panel shows the distribution of "luminous" satellites (see § 2.2), and the right panel 
shows the distribution for the 11 subhalos with the largest V max in each host. In each panel, the dotted lines represent halo Gi, the dashed lines represent halo G2, 
and the dot-dashed lines represent halo G 3 . The thin, solid lines in each panel represent the expected CDF for an isotropic satellite distribution. 



and find best-fit planes to the satellite positions in our host ha- 
los. Specifically, we determine a best-fit plane by minimizing 
the root-mean-square (rms) of the perpendicular distances of all 
satellites to the plane, 

D -=V^h/ — • (2) 

In Eq. (2), h is a unit vector normal to the plane, d is the per- 
pendicular distance from the origin (in this case, the origin is 
the Galaxy or the halo center) to the plane, x^ is the position of 
the ^ith satellite, and the sum is over all N satellites. We define 
a polar angle < u> < it/ 2, as the angle between a satellite po- 
sition vector and the best-fit unit vector n, set at the point on 
the plane that minimizes the distance to the origin. Addition- 
ally, the value of D rms itself can also be used as a measure of 
the planarity of the satellite distribution (Kang et al. 2005), and 
we consider this statistic below. 

3. RESULTS 
3.1. Visual Impression 

Figure 2 gives a visual impression of the planarity of the 
observed MW satellites, as well as the subhalos and the lumi- 
nous subhalos in our simulations. Here and below, we consider 
all eleven known MW satellites within R max = 300 kpc (Grebel 
et al. 2003). The MW satellites are shown in the rightmost pan- 
els of Fig. 2. The planarity of the distribution, noted by Kroupa 
et al. (2005) is very clear. Note, however, that recent studies 
indicate that the poles of the orbits of Ursa Minor (Piatek et al. 
2005), as well as Sculptor and Fornax (Piatek et al., in prepa- 
ration), are not coincident with the pole of the plane of MW 
satellites, nor are they coincident with each other. This is an 
indication that the high degree of planarity of the MW dwarfs 
may be transitory. 

In the top row of Fig. 2, we show projections of all subhalos 
of each host onto the plane perpendicular to their best-fit plane. 
In these projections the best-fit plane runs vertically along the 
axis x = and we show the projected direction of the major 
axes in each panel by a thin solid line. In the bottom row, we 



show similar projections for the luminous subhalos. The spatial 
distributions of the highest- V max subhalos are similar to those of 
the luminous subhalos. 

The most obvious feature in all of the panels is that the best- 
fit plane of subhalos is strongly correlated with the major axis 
of the DM distribution of the host halo. If we draw 10 5 random 
permutations of eleven subhalos from the full sample, the prob- 
ability that the best-fit plane of any such subsample is inclined 
by more than 45° from the host halo major axis is only P ~ 9%, 
signaling anisotropy of the subhalo population as we discuss in 
more detail below. Correspondingly, the best-fit planes for the 
luminous subhalos are also aligned with the major axes of their 
hosts. Their degree of planarity is visually comparable to that 
of the MW satellites. In the following subsections, we quantify 
the anisotropy of satellite distribution. 

3.2. Satellite Planarity 

Kroupa et al. (2005) argued that the planarity of the observed 
MW satellites is inconsistent with the CDM paradigm based on 
their distribution in |cos(w)|. However, their conclusion is in- 
correct for two reasons. First, as we show below, CDM does 
not predict that dwarf-sized subhalos are distributed isotropi- 
cally. Second, the null hypothesis for the isotropic distribution 
for I cos(w)| used by Kroupa et al. (2005) is correct only in the 
limit of large sample size. This is because for a small num- 
ber of objects, the distribution of |cos(w)| is not related to a 
fixed reference frame, but a frame determined by the selected 
objects themselves. This modifies the underlying cumulative 
distribution function (CDF) for |cos(cj)|. Kroupa et al. (2005) 
computed the distribution of ui for a single sample of 10 5 sub- 
halos, the result of which is a nearly uniform distribution in 
|cos(w)| on the interval [0, 1]. However, the |cos(o;)| CDF for 
small number of objects drawn from the isotropic distribution 
is different. For example, consider the limit of randomly se- 
lecting precisely three objects from any underlying distribution. 
The three objects will always lie on the best-fit plane, such that 
I cos(w)| = for all objects, independent of their underlying dis- 
tribution. 

In order to account for the small sample size, we generate the 
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FIG. 4. — The cumulative fraction of satellites with the absolute value of the cosine of the zenith angle < | cos(#)|. The zenith angle, < 8 < n, is defined as the 
angle from the major axis of the DM distribution of the host. The left panel shows the distribution for all subhalos with V max > 12kms~' . The center panel shows 
the distribution of only those subhalos that are "luminous" according to the model of Kravtsov et al. (2004, see § 2.2). The right panel shows the distribution for the 
eleven subhalos with the largest V max in each host. The observed MW satellites (thick solid line) are placed on this plot by defining 6 to be the angle from the pole 
of the MW disk. In each panel, the dotted lines represent halo G i , the dashed lines represent halo G2 , and the dot-dashed lines represent halo G3 . 



distribution in | cos(w) | that should be expected from drawing 
eleven satellites from a uniform distribution. First, we assume 
an underlying isotropic subhalo distribution with a radial num- 
ber density profile n(r) oc 1 /(l +(r/r c ) 3 ) with r c = 0.25/?i8o (see 
Zentner et al. 2005, who found this to be a good description 
of subhalo distributions in simulations and their semi-analytic 
models). Our results are not sensitive to the particular as- 
sumptions that we make about the radial distribution. We then 
draw 10 5 random samples of eleven satellites from this distri- 
bution, determine the best-fit plane for each sample, and com- 
pute cos(w) for each satellite and the corresponding CDF. The 
resulting average CDF of |cos(dj)| is shown as the thin, solid 
line in Figure 3. Notice that the | cos(w)| distribution for eleven 
satellites selected from an isotropic distribution is significantly 
different from a uniform distribution in |cos(w)|, which would 
be a diagonal line. 

We can assess the probability that the observed MW satel- 
lites are drawn from an isotropic distribution using the CDF 
of |cos(w)| constructed as described above. The Kolmogorov- 
Smirnov (KS) probability to draw the observed distribution of 
MW satellites within R < 300 kpc from an isotropic distribution 
isP^ s s °~0.15. 6 

In a similar analysis, we draw 10 5 random samples of eleven 
subhalos from the full, simulated subhalo populations (V max > 
12kms _1 and R < 300 kpc) and construct the CDFs for these 
samples as shown in Fig. 3. Comparing this distribution to the 
isotropic case gives a KS probability of P K s — 0.49. However, 
as we noted above, the best-fit planes in the subhalo samples are 
inclined by less than 45° with respect to the major axis of the 
host halo in 91% of the samples, while there is no such prefer- 
ential alignment in the isotropic case. For the luminous subha- 
los, the KS probability that they are drawn from an isotropic 
distribution is Pks — 0.39. Finally, the probability that the 
MW satellites are drawn from the distribution of all subhalos is 
P K s — 0.16. Figure 2 and the analysis below show that the dis- 

6 If we were to follow Kroupa et al. (2005) and use the uniform CDF that results 
from a very large sample size, we find P^® 2 X 10~ 3 , in good agreement with 
their analysis. 



tributions of both the MW satellites and subhalos are actually 
anisotropic. The conclusion is that for a small sample (eleven), 
the |cos(o;)| distribution is non-discriminatory. Other statistics 
must be used to draw meaningful conclusions. 

The problem of the reference frame uncertainty for small 
sample sizes can be avoided if the frame is fixed and does not 
depend on the objects. In particular, we can choose to define 
the angles with respect to the halo major axes, (i.e., the zenith 
angle 9, defined in § 2.3). We show the CDF of |cos(6>)| for 
the simulated subhalos in the left panel of Figure 4 and the 
CDF for an isotropic distribution by the thin diagonal line. 
For all three of the simulated MW-sized hosts, the subhalos 
are clearly distributed anisotropically with a preferential align- 
ment along the host halo major axis. The KS probability of 
drawing this subhalo population from an isotropic distribution 
is P^ s s ° ~1.5x 10~ 4 . This is in contrast to the result of Willman 
et al. (2004), who find an isotropic distribution of subhalos in 
a MW-sized halo simulated by Reed et al. (2003). We find the 
differential fraction of satellites per unit x = | cos(#)| averaged 
over all three host halos to be sharply peaked and well described 
by 

^^0.80+1.15x 475 . (3) 
dx 

The thick, solid lines in each of the panels of Figure 4 show 
a similar distribution in zenith angle for the eleven observed 
MW satellites within 300 kpc. The orientation of the MW halo 
is unknown, so the MW satellites cannot be placed on this plot 
without making an assumption about the orientation of the halo 
with respect to the disk. In the reference frame in which the 
zenith angle 9, is measured with respect to the Galactic pole 
(we discuss implications of this assumption in § 4), the KS 
probability of drawing the MW satellites from a distribution 
that is isotropic about the Galactic center is P^ s s ° ~ 0.07, while 
the probability of drawing the MW satellites from the sample 
of CDM subhalos is Pks — 0.34. Alternatively, if we assume 
that the MW satellites are aligned along the minor or interme- 
diate axes, the probabilities for the subhalos and MW satel- 
lites to have the same parent distribution are P K s — 0.02 and 
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Pks — 0.03, respectively. 

For the subhalos that are luminous according to the KGK04 
model, the preferential distribution at low zenith angles re- 
mains and the KS probability of for all luminous subhalos 
to be drawn from an isotropic underlying distribution is only 
P^ s s ° ~ 0.02. For the eleven highest-V max subhalos in each host, 
the KS probability of being sampled from isotropic distribution 
is P^ s s ° ~ 0.05. Note that for such small samples of objects, the 
halo-to-halo variation in the zenith angle distribution is large. 

A useful measure of planarity is the rms distance to the best- 
fit plane, D rms (Kang et al. 2005, see Eq. (2)). In agreement 
with Kroupa et al. (2005), we find that ~ 26.4kpc for 
the eleven MW satellites within 300kpc. For eleven satellites 
selected randomly from the isotropic distribution above, the 
mean value of the rms plane width in 10 5 sample realizations 
is Df° s ~ 72kpc and the probability of a distribution similar to 
the MW satellites is P(D™ S < D^f) ~ 5 x 10" 3 . Similarly, for 
the subhalo realizations averaged over all three host halos, the 
mean width is D lms ~ 58kpc and P(D rms < D^) ~ 0.02. This 
shows that both MW satellites and DM subhalos are distributed 
anisotropically, but the observed distribution is somewhat more 
planar. 

The radial distribution of observed satellites is, however, 
more centrally concentrated than that of the subhalos. The me- 
dian distance of the MW satellites from the Galactic center is 
a factor of two smaller than the median distance of subhalos 
from the centers of their hosts (Taylor et al. 2003; Kravtsov 
et al. 2004). The luminous subhalos have a radial distribution 
that is similar to the MW satellites and represent a biased sub- 
sample of the overall subhalo population. The rms distance to 
the best-fit plane is sensitive to the radial extent of the popula- 
tion and we expect more centrally concentrated populations of 
objects to have smaller values of £> rms . 

Indeed, we find that the peak of the distribution of D rms scales 
linearly with the median of the radial positions of the underly- 
ing population /? me d- The rescaled rms distance to the plane, 
5 ms = Z) rms //? me d, can serve as a useful measure of planarity 
that normalizes out the radial extents of the different popula- 
tions. In Figure 5, we show the probability distribution of se- 
lecting eleven satellites with a particular value of Srms drawn 
from the isotropic distribution and from the samples of subha- 
los of hosts Gi, G2, and G3. For both the isotropic distribu- 
tion and for the simulated subhalos, /? me( j « 167kpc. Figure 5 
shows that in terms of the probability distribution of <5 rms , the 
observed value is not unusual. In fact, we find that the value of 
D^f is not unlikely (P(<5^ s < 6%™) ~ 0. 14) even for a random 
sample of eleven objects drawn from an isotropic distribution, 
if their radii are rescaled to have the same R me d as the MW 
satellites. The figure also shows that the distribution of sub- 
halos is anisotropic, which shifts the probability distribution of 
Jrms, increasing the probability of selecting a subsample with 
a degree of planarity comparable to that of the MW satellites: 
P(<Cs<O-0-30. 

In the model of KGK04, the luminous subhalos of Gi and 
G3 match the observed radial distribution of the MW satellites 
well. We thus expect that their D rms will be comparable to the 
observed value. Indeed, we find for the luminous subhalos that 
D rms ~ 29.4kpc, D rms ~ 46.2kpc, and D lms ~ 27.4kpc for hosts 
Gi, G2, and G3, respectively. These values are similar to Djmf . 
The values of <5 rms for the MW satellites and the luminous sub- 
halos are shown by arrows in Figure 5. 

These results indicate that the main reason that the MWsatel- 
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FIG. 5. — The distribution of the rms dispersions of satellite distribution 
around their best fit planes. The dispersions are rescaled to the median distance 
of each population to the center of their host (see text for details). The dotted 
line corresponds samples of 11 objects selected from an isotropic population 
distributed radially as the simulated subhalos. The other lines corresponds to 
samples of 1 1 subhalos drawn from the subhalos of Gi (solid line), G2 (dashed 
line), and the G3 (dot-dashed line). The labeled arrows at the bottom of the 
plot mark the values of Drms/Rmed' f° r me observed MW satellites and the 
luminous subhalos of Gi, G2, and G3. 



lites occupy a narrower plane than DM subhalos is their more 
centrally-concentrated radial distribution. Given that there are 
physical reasons to expect such radial bias (KGK04), we con- 
clude that the observed anisotropy of the MW satellites is con- 
sistent with CDM predictions. However, we note that this 
agreement requires approximate alignment of the major axis 
of the halo that hosts the MW and the pole of the MW disk. In 
the next section, we discuss the origin of the subhalo anisotropy 
and the implications of such disk-halo alignment. 

4. DISCUSSION 

4.1. Origin of the Subhalo Anisotropy 

It is interesting to ask why satellite halos should have a 
strongly-anisotropic distribution aligned with the major axis of 
the host halo. 

Several suggestions have been made to explain the 
anisotropy observed by Holmberg (1969) and Zaritsky et al. 
(1997). Quinn & Goodman (1986) investigated the effect of 
enhanced dynamical friction for orbits that are nearly co-planar 
with a galactic disk. The idea is that nearly co-planar orbits 
would be driven toward the disk plane (e.g., Binney 1977) and 
decay more rapidly due to additional interactions with the disk 
component. In this way, satellites on nearly co-planar orbits 
would be preferentially cannibalized by the disk and these or- 
bits depopulated. For example, such cannibalized dwarf galax- 
ies can significantly contribute to the formation of the thick disk 
(Abadi et al. 2003). However, the conclusion of Quinn & Good- 
man (1986) was that this process is not efficient enough to ac- 
count for the results of Holmberg (1969) for satellites closer 
than ~ 50kpc. Penarrubia et al. (2002) extended this argument 
to include the effect of an oblate DM halo. They found this 
process to be efficient only within ~ 50kpc of the host due to 
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FIG. 6. — Trajectories of 100 randomly selected subhalos accreting onto host halos Gi, G2, and G3 (panels from left to right). The trajectories are plotted relative to 
a coordinate system centered on the most bound particle in the parent halo and oriented along the principal axes of the parent halo evaluated at z = 0. The coordinates 
are comoving. The first point of each trajectory is at the time when the comoving distance between subhalo and host first became smaller than 600/i kpc. In each 
panel, the ellipsoids have major axes of length f = 290/r'kpc, and the axis ratios are set to the axis ratios of each host at £ = Rigg and z = 0. The thick, solid, 
mutually-orthogonal lines denote the principal axes of the host halos. We label each coordinate axis by the corresponding principal axis (a> b > c) with which it is 
aligned. Preferential accretion along the direction of the long axis is evident for halos Gi and G3. 



considerably longer orbital decay times at larger radii (Zarit- 
sky & White 1994; Zaritsky & Gonzalez 1999). Note that the 
MW-sized halos in our simulations are prolate, so this process 
should be even less efficient. Yet, the anisotropy of satellites is 
present. 

Two possibilities are that the anisotropy reflects a direction of 
preferential infall set by the environment (e.g., Tormen 1997) 
and/or that there is some other dynamical process that drives 
evolution toward anisotropy after accretion onto the host halo. 
Knebe et al. (2004) investigated the first possibility for clus- 
ter halos to address the anisotropy observed in systems such as 
the Virgo cluster (e.g. West & Blakeslee 2000). Knebe et al. 
(2004) concluded that preferential accretion along the direc- 
tions of filaments accounts for much of the bias in satellite or- 
bits in cluster-sized systems. This result does not extend triv- 
ially to MW-sized systems. Clusters are rare, highly-biased ob- 
jects that generally form relatively recently at the "nodes" of 
filaments that are comparable in thickness to the size of clus- 
ters themselves (e.g., Klypin & Shandarin 1983; Bond et al. 
1996; Colberg et al. 2004). It is therefore not surprising that 
many merging halos get "funneled" into clusters along domi- 
nant filamentary directions. Unlike clusters, MW-sized objects 
typically form earlier, are significantly less biased at the present 
epoch, and do not generally reside in filaments of comparable 
dimension. In this case, it seems that infalling substructure may 
be less likely to accrete along a single, dominant filamentary di- 
rection over the formation history of the halo. In addition, clus- 
ters have generally assembled their masses and accreted their 
satellites much more recently than galaxy-sized systems. As a 
result, many cluster subhalos are dynamically young and have 
undergone only a small number of orbits (< 1-2) in the poten- 
tial of the main halo, while the subhalos in our MW size halos 
have typically undergone several orbits within the main halo 
(see, e.g., KGK04). 

In Figure 6, we show trajectories for 100 randomly-selected 
subhalos of each host halo as they are accreted. The trajecto- 
ries are constructed from the 96 saved simulation timesteps by 



examining the 25% most bound particles in each halo. Halos 
at adjacent timesteps that share the highest common fraction of 
these particles are identified with each other, as described in § 4 
of KGK04. In two of the three cases (halos Gi and G3), there 
is clearly preferred satellite accretion along the direction of the 
major axis of the host halo. In the case of halo G2, the accretion 
is less anisotropic, but there appears to be a small preference for 
accretion toward the octant a < 0, b < 0, c> 0. We have exam- 
ined the accretion histories of the host halos more closely and 
halo G2 appears to feed off of a filament that runs roughly in 
this direction during its early evolution (z> 1). It subsequently 
accretes along a filament that is more closely aligned with its 
major axis direction at z = 0. Note that the halos Gi and G3 
have both the highest degree of satellite anisotropy (see § 3) 
and the most pronounced preferred accretion directions. 

These results indicate that a preferential direction of satellite 
accretion is an important factor in determining the distribution 
of satellites in MW-sized halos, just as it is for cluster-sized 
hosts. Although filaments at the present time may be thick, they 
were significantly thinner in the past, when many of the surviv- 
ing subhalos were accreted. Moreover, the matter distribution 
in filaments is concentrated toward the axis of the filament (e.g., 
Colberg et al. 2004). Finally, a fraction of subhalos may be ac- 
creted as members of groups, which are biased spatially and 
often located near the centers of filaments. That the preferen- 
tial accretion direction is correlated with the major axis of the 
host halo is not surprising, because the major axis is typically 
determined by the direction of the most recent major merger. 

In our DM-only simulations, effects like the cannibalization 
of satellites on co-planar orbits by material associated with the 
disk are absent, yet the substructure anisotropy is present, even 
for subhalos that have orbited within the host potential for many 
dynamical times. This is due to the fact that elongated poten- 
tials, similar to the potentials induced by our triaxial host halos, 
support orbits that make long excursions along the major axis 
of the potential (see Statler 1987). To illustrate this, we perform 
a simple, idealized experiment. We integrate the orbits of 200 
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FIG. 7. — Generation of an anisotropic distribution of satellites in triaxial potentials. The distribution of zenith angle cosines for 200 isotropically accreted test 
particles at different epochs (lines of various type; the time is in units of the dynamical time of the host at the virial radius). The left panel represents control case 
of the spherical halo potential, while in the right panel the host has ellipsoidal isodensity contours with axis ratios of b/a = 0.65 and c/a = 0.60, similar to those of 
halos Gi, G2, and G3. The fiducial host halo has a virial mass of M 180 = 10 12 r'M and the NFW mass distribution with concentration parameter of ciso = 12. 
The solid line indicates an integration for 1 ti yn , the dot-short-dashed line indicates integration for 3f(i yn , the dashed line for 5(dyn, the dotted line for 15<j yn , and 
the dot-long-dashed line for 50fj yn . The thin solid line corresponds to the isotropic distribution. The figure shows that in the spherical potential tracer distribution 
remains isotropic at all times. In this case of prolate, triaxial potential, the orbits evolve toward an anisotropic distribution similar to that of the observed satellites 
and simulated subhalos. 



test particles in the static potential of a triaxial generalization of 
the density profile of Navarro et al. (1997, NFW hereafter), 

' ,oo< (£)'( i+ £) 2 ' <4) 

with ellipsoidal contours of constant density. ( is defined as in 
§ 2.3. For simplicity, we choose constant axis ratios b/a = 0.65 
and c/a = 0.6, similar to our simulated host halos (Fig. 1), and a 
halo concentration c\%q = Ri%o/r s = 12, which is typical of MW- 
sized halos (Bullock et al. 2001). We choose the initial veloci- 
ties for each orbit according to the distribution of initial condi- 
tions for subhalos presented in Zentner et al. (2005). In order to 
demonstrate the influence of the triaxial potential, we assumed 
spherically symmetric infall, rather than the anisotropic infall 
depicted in Figure 6. 

The results of this exercise are shown in Figure 7, where we 
plot the distribution of subhalo zenith angle for the subhalos in 
the triaxial model and a spherical model that serves as a con- 
trol. We integrate the orbits over an interval of 50fd yn , where 
fdyn is the dynamical time of the model at R\%q. The net effect is 
clear and not surprising, but is often neglected. The prevalence 
of orbits that make long excursions along the major axis of the 
potential induces and maintains an anisotropic distribution of 
test particles. The test particles assume a distribution consistent 
with the triaxiality of the potential and, in fact, in a live, triaxial 
density distribution, it is the presence of these orbits that main- 
tain the triaxial shape of the system (Gerhard & Binney 1985; 
Udry & Martinet 1994; Barnes & Hernquist 1996; Merritt & 
Quinlan 1998; Valluri & Merritt 1998). 

We neglected many effects in this experiment. These in- 
clude the influence of the subhalos on the host potential, and 
the growth and evolution of the host halo potential that occurs 
in cosmological simulations. As we show in Figure 1, the host 



halo shapes evolve from c/a^ 0.4 at 1 to c/a^0.6 atz^O. 
Moreover, subhalos on elongated orbits will be more vulnerable 
to tidal disruption because they generally pass closer to the cen- 
ter of the potential where the tides are strongest (e.g., Zentner 
& Bullock 2003; Zentner et al. 2005). The net effect should be 
that such orbits will be gradually depopulated. Nevertheless, 
our experiment illustrates that regardless of anisotropic infall, 
the satellite distribution can develop anisotropy due to the triax- 
iality of the host potential, even in the absence of any enhanced 
destruction due to the central disk. 



4.2. Anisotropy of the globular cluster distribution 

As we pointed out above, our results indicate that the ob- 
served distribution of the MW dwarf satellite galaxies is gen- 
erally consistent with CDM predictions if the major axis of the 
MW DM halo is approximately aligned with the Galactic pole. 
If this is the case, based on our discussion of orbital structure in 
the previous section, we should expect similar anisotropics to 
exist in other populations that may serve as test particles in the 
MW halo potential. 

Consider the distribution of the MW globular clusters (GCs). 
Frenk & White (1982) studied a subsample of the MW GCs and 
concluded that both metal-rich and metal-poor globular cluster 
systems are slightly flattened. Subsequently, Zinn (1985) used 
a larger GC sample and concluded that metal-rich GCs are in 
a disk-like configuration (see also Armandroff 1989), while the 
distribution of the metal-poor clusters is nearly spherical. More 
recently, Hartwick (2000) analyzed the distribution of 15 metal- 
poor, distant globular clusters (/?> 25kpc and [Fe/H] < -1) 
and found that they form a flattened system with a minor axis 
highly-inclined relative to the MW disk rotation axis. 

Following Zinn (1985), we divide the GC sample of Harris 
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FIG. 8. — The distribution of metal-poor, MW globular clusters ([Fe/H] < 
-0.8) and dwarf satellite galaxies. All positions are projected onto a plane per- 
pendicular to the best-fit plane of the 44 most distant globulars (R > 10 kpc). 
The points correspond to the GCs with R < lOkpc. The open triangles corre- 
spond to GCs with R > 10 kpc. For comparison, we show the positions of the 
MW satellite galaxies as the filled squares. The thin, solid, horizontal line rep- 
resents the orientation of the MW disk. Note that only 8 of the MW satellites 
fall within the limits of this plot and the projection is oriented at an angle of 
nearly 36° with respect to their best-fit plane. 



(1996) 7 into halo and disk clusters using a metallicity threshold. 
We consider the distribution of metal-poor ([Fe/H] < -0.8), 
halo GCs. The distribution of the metal-poor GCs is shown in 
Figure 8 along with the distribution of the innermost MW dwarf 
satellites. Visually, the distribution of these two populations of 
halo objects appear to have similar anisotropy. 

Just as for the dwarf satellites, we compare the distribution 
of GCs to an isotropic distribution using the distribution of 
|cos(f9)|, where 8 is defined relative to the Galactic pole, and 
<5^ because |cos(w)| has little discriminatory power. The re- 
sults are summarized in Table 4.2. The inner globular clusters 
(R < 10 kpc) appear to lie in a flattened distribution (P(6™ s < 
6^ s ) ~ 0.02) that is closely aligned with the disk plane. The 
normal to the best-fit plane of inner halo GCs is offset form the 
MW pole by only 6 M w ^17°. Similar to Hartwick (2000), 
we find that the distant halo GCs show marginal evidence for 
a distribution that is aligned with the pole of the MW disk 
(©MW — 80°), a configuration that is consistent with the dis- 
tribution of the MW dwarf satellites. 

An additional piece of evidence comes from the distribution 
of the M31 satellites. In this case, distances to the M31 sub- 
group have considerably larger errors than the MW dwarfs, 
so these objects carry less statistical weight. Nevertheless, 
Hartwick (2000) studied the satellites of M31 in detail and 
found that, like the MW satellites, they appear to be arranged 
in a prolate structure that has its best-fit major axis misaligned 
with the pole of the M31 disk by only ~ 20°. The distributions 
of M3 1 and MW satellites appear to show a similar anisotropy 
with preferential alignment along the galactic poles. 

4.3. Implications 

7 The globular cluster catalog is available at URL 

http : //physwww . physics . mcmaster . ca/~harris/mwgc . dat. 
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Column description: (1) the range of galacto-centric distances in GC 
subsample; (2) Omw, the angle between the MW pole and the unit 
vector normal to the best-fit plane of GCs in the sample (6mw = deg 
corresponds to the Galactic disk plane, while 6mw = 90 deg corre- 
sponds to a plane parallel to the MW pole); (3) Gsats, the angle be- 
tween the best-fit plane of GCs and the best-fit plane of MW dwarf 
satellite galaxies; (4) P| s , the KS probability for the GC sample to be 
drawn from an isotropic distribution using the cumulative distribution 
of |cos((9)|, where 9 is the angle with between position vector of the 
satellite and Galactic pole; (5) the number of halo GCs in each sub- 
sample. 

The angular distribution of satellites and its interpretation in 
the context of the anisotropy of the CDM subhalos have im- 
portant implications for our understanding of galaxy forma- 
tion. Collapsing halos acquire angular momentum via interac- 
tions with quadrupolarmass density fluctuations (Peebles 1969; 
Doroshkevich 1970; Efstathiou & Jones 1979; Barnes & Efs- 
tathiou 1987). In the simplest scenario of disk galaxy forma- 
tion, baryons in halos begin by sharing the angular momentum 
distribution of the DM, on average, conserving it as they cool 
and condense (e.g., Fall & Efstathiou 1980). This leads to a pic- 
ture where the poles of disk galaxies are collinear with the net 
angular momentum vectors of their host halos, which are gen- 
erally aligned with halo minor axes (e.g., Warren et al. 1992; 
Porciani et al. 2002; Faltenbacher et al. 2005). 

van den Bosch et al. (2002) studied the alignment of the an- 
gular momenta of DM halos and their baryons in adiabatic sim- 
ulations, and found that they are generally well aligned, with a 
distribution for the angular misalignment that is sharply peaked 
between ~ 10° -20°, but with an extended tail to larger angles 
such that the median is ~ 30°. Chen et al. (2003) found similar 
results in adiabatic simulations and models with radiative cool- 
ing. Both van den Bosch et al. (2002) and Chen et al. (2003) 
measured increasing misalignment with decreasing halo mass 
and speculated that the extended tail could be partially due to 
noise in the measurement of angular momenta in small objects. 
We note that Kazantzidis et al. (2004a) found the gaseous disk 
of a young (z> 4) galaxy progenitor in their simulation to be 
aligned nearly perpendicular to the major axis of the DM halo, 
which is the alignment that our results suggest. 

In this paper, we have shown that the predicted spatial dis- 
tribution of CDM subhalos is consistent with the distribution 
of the MW satellites, if the pole of the MW disk is nearly 
aligned with the major axis of its outer DM halo. This requires 
a more complicated disk formation scenario where the halo and 
baryons mutually adjust as they evolve toward a stable configu- 
ration. It would be interesting to explore whether such an align- 
ment is supported by observational evidence from other galax- 
ies and numerical models of galaxy formation. Along these 
lines, we show the anisotropic, two-dimensional projected sub- 
halo distributions in the Appendix. 
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As we noted in the introduction, current observational results 
are contradictory. Holmberg (1969), Zaritsky et al. (1997), and 
Sales & Lambas (2004) report that satellites of other galaxies 
exhibit statistical anisotropy similar to that of the MW dwarfs. 
Specifically, they find satellites to be preferentially located near 
the minor axes of the projected light distributions of host galax- 
ies. However, the study of Brainerd (2004) shows evidence of 
the opposite correlation of satellite position with the major axis 
of their host galaxy. 

One potential avenue for checking the disk-halo alignment 
is weak lensing. Hoekstra et al. (2004) recently presented the 
first weak lensing measurement of halo ellipticity obtained un- 
der the assumption that the halo mass and galaxy light distribu- 
tions are aligned. This is likely a sound assumption for early- 
type galaxies, which may indeed dominate the lensing signal in 
a sample of mixed morphological types. However, if galactic 
disks are preferentially aligned orthogonal to the major axes of 
their halos, the alignment for late-type systems would be oppo- 
site to that assumed by Hoekstra et al. (2004). Observational 
tests should be possible with the large dataset of the Sloan Dig- 
ital Sky Survey. 

On the theoretical side, Navarro, Abadi, & Steinmetz (2004) 
recently argued that the angular momenta of galactic disks are 
aligned perpendicular to the minor axis and parallel to the in- 
termediate axis of the inertia tensor of the surrounding matter 
distribution, both in their simulations of galaxy formation and 
in observations of nearby galaxies in the Local Supercluster. 
Unfortunately, it is not clear whether the major axes of galac- 
tic halos are always aligned with the intermediate axis of the 
surrounding large-scale structure. Navarro et al. (2004) also 
reported that the disks of nearby spirals tend to be oriented per- 
pendicular to the Supergalactic plane. We generally expect the 
major axes of halos to be oriented along the filamentary struc- 
ture of the Local Supercluster, so such a tendency is consistent 
with a scenario where disks orient themselves perpendicular to 
halo major axes. 

The shape of the halo and the orbits of satellites within pro- 
late and triaxial hosts may have important consequences for 
disk heating and the build up of the stellar halo. For exam- 
ple, given that the orbital energies of most surviving subhalos 
are relatively large, one may expect that the number of satellite 
passages close to the galactic center within the last ~ lOGyr 
should be small (e.g., Font et al. 2001). However, a satellite on 
a box-like orbit can pass arbitrarily close to the disk. This may 
enhance both the heating of the stellar disk (e.g., Toth & Os- 
triker 1992; Huang & Carlberg 1997; Ibata & Razoumov 1998; 
Moore et al. 1999; Velazquez & White 1999) and the efficiency 
of tidal disruption of satellites and the formation of the stellar 
halo (e.g., Helmi & White 1999; Bullock et al. 2001). Objects 
that experience close passages will be preferentially disrupted 
during halo evolution, such that it may appear that the present- 
day halo is devoid of satellites capable of doing damage. These 
considerations imply that the triaxiality of host potentials has 
to be taken into account in semi-analytic and numerical models 
of subhalo evolution (e.g., Bullock et al. 2000; Taylor & Babul 
2001, 2004; Benson et al. 2004; Zentner & Bullock 2003; Zent- 
ner et al. 2005) in order for the treatment to be accurate. 

Lastly, the alignment of satellite orbital planes along the 
halo major axis can have implications for the evolution of tidal 
tails. Helmi (2004a,b) argues that the data on the Sagittarius 
tidal debris (Majewski et al. 2004; Law et al. 2004) provide a 
strong indication that the MW halo is prolate with its major 



axis aligned with the Galactic pole and that the mean axis ratio 
within the orbit of Sagittarius is 0.65 < c/a< 0.8. Again, this 
is the alignment necessary to explain the positions of the MW 
dwarfs. Note, however, that Johnston et al. (2004) argue that 
the precession of the orbital plane of Sagittarius implied by the 
data require an oblate halo that is, at most, slightly flattened 
c/a> 0.85. 

4.4. Caveats 

One of the caveats to our results (and, indeed, to the results 
of the related studies of Knebe et al. (2004) and Kang et al. 
(2005)) is that the simulations that we present follow the dissi- 
pationless evolution of DM only. The net effect of dissipation 
and the condensation of baryons on the anisotropy of satellites 
is unclear. There are several aspects to consider. 

Cosmological gasdynamics simulations of galaxy and cluster 
formation show that radiative gas cooling results in DM halos 
that are significantly more spherical in their inner regions than 
halos in dissipationless simulations (Kazantzidis et al. 2004a). 
At Q = 0.1^i8o, the average increase in the minor-to-major axis 
ratio is A(c/a) ~ 0.3, but this shift is a declining function of 
C, such that at C> 0.5 it is A(c/a)< 0.1. This does not affect 
the accretion of satellites along preferred directions (Fig. 6), but 
will make the halo potential in the inner regions more spheri- 
cal, reducing any alignment of orbits with the halo major axis. 
However, this effect is small at large galacto-centric radii and 
satellites with large apocenters may still move under the influ- 
ence of an effectively prolate potential. 

The presence of a disk could cause the DM halo in its vicin- 
ity to adopt a locally-oblate shape aligned with the disk plane 
(as opposed to the generally prolate shapes seen in cosmologi- 
cal A^-body simulations), again enhancing destruction of satel- 
lites on prograde, co-planar orbits due to anisotropic dynamical 
friction. If the disk is oriented perpendicular to the halo major 
axis this would enhance satellite anisotropy compared to our re- 
sults. Finally, it is possible that the observed, strong anisotropy 
of the MW satellites is partially due to obscuration by the disk. 
The incompleteness of the current satellite sample is uncertain, 
though Willman et al. (2004) argue that a significant fraction of 
MW satellites at large distances may still be undetected. 

5. SUMMARY AND CONCLUSIONS 

We study the spatial distribution of dwarf subhalos in MW- 
sized DM halos using dissipationless cosmological simulations 
of the concordance flat ACDM cosmology. Specifically, we 
compare the simulated subhalo populations with the observed 
distribution of the known MW satellite galaxies. We also test 
whether the predictions of CDM simulations are consistent with 
observations using two possible scenarios for mapping lumi- 
nous satellite galaxies onto subhalos in MW-sized DM halos. 
In the first scenario, the luminous dwarf satellites are identified 
using the semi-analytic model for dwarf galaxy formation pro- 
posed by Kravtsov et al. (2004). In the second, the luminous 
satellites reside in those few subhalos with the largest maxi- 
mum circular velocities V max (Stoehr et al. 2002, 2003). Our 
main results can be summarized as follows. 

1. The distribution of subhalos in host DM halos is not 
isotropic. Subhalos are preferentially aligned with the ma- 
jor axis of the triaxial host halo mass distribution. The 
KS probability of choosing the subhalo populations in our 
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simulated MW-sized host halos from an isotropic distri- 
bution is 1.5x10^. 

2. The method used by Kroupa et al. (2005) does not demon- 
strate that the MW satellites are inconsistent with either an 
isotropic underlying distribution or a distribution similar 
to that of subhalos in CDM simulations. We argue that 
Kroupa et al. (2005) adopted an incorrect null hypothesis. 

3. In terms of angular distribution, the probability for the 
MW satellites to be drawn from the anisotropic distribu- 
tion of simulated subhalos is Pks — 34% under the as- 
sumption that the MW pole is aligned with the major axis 
of its host halo. Alternatively, if the MW pole is aligned 
with the minor axis of the host halo, the KS probability is 
Pks ^ 2%. 

4. The apparent planarity of the MW satellite distribution 
can be explained by the anisotropy of the subhalo dis- 
tribution and the relative radial bias of luminous dwarf 
satellites relative to subhalos. Specifically, the subhalo 
subsamples that correspond to the luminous subhalos in 
the model of KGK04 and the eleven highest- V max subha- 
los exhibit a degree of planarity that is similar to the ob- 
served MW satellites. Note again that in our simulations 
such planar distributions are likely to be nearly aligned the 
major axis of the host halo. This, in turn, implies that near 
alignment of disk pole and halo major axis is required to 
explain the observed satellite distribution. 

5. In agreement with Hartwick (2000), we find that distant 
(R> 10 kpc), metal-poor ([Fe/H] < -0.8), MW globu- 
lar clusters exhibit anisotropy similar to that of the dwarf 
satellites. 

6. The observed anisotropy of the MW satellites compared 
with the CDM predictions for subhalo orientations, along 
with evidence for the Holmberg effect in other galaxies 
(Holmberg 1969; Zaritsky et al. 1997; Sales & Lambas 
2004) including the dwarf satellites of M31 (Hartwick 
2000), the distribution of MW halo globular clusters, and 
the indirect arguments of Helmi (2004b), provide evi- 
dence for a consistent picture in which the outer DM halos 
surrounding spiral galaxies should be nearly perpendicu- 
lar to the disk planes. This has interesting implications for 
the understanding of disk galaxy formation and the orbital 
evolution of satellite galaxies. 

We discuss the origin of the anisotropy of the subhalo distri- 
bution in simulations and show that, similar to galaxy clusters, 
Galaxy-sized halos accrete substructure along a preferential di- 
rection that is strongly correlated with the major axis of the host 
halo (Fig. 6). We also stress that orbital evolution in a triaxial 
potential results in an anisotropic spatial distribution of tracer 
objects, even if their accretion is isotropic. 

The fact that consistency of observations with the ACDM 
prediction appears to require near alignment of disk angular 
momenta with the major axis of host halos is surprising. The 
angular momenta of halos are typically aligned with their minor 
axes (e.g., Warren et al. 1992; Porciani et al. 2002; Faltenbacher 
et al. 2005). This could indicate that disk formation is accompa- 
nied by the evolution of the angular momentum of the baryonic 
material as it cools and condenses in the center. This would not 



be surprising given that condensation of baryons has a signif- 
icant impact on the density and shape of the surrounding DM 
halo (Gnedin et al. 2004; Kazantzidis et al. 2004a). However, 
more detailed numerical studies of disk galaxy formation are 
needed to understand the underlying processes in detail. 

It would be interesting to test whether the putative disk-halo 
configuration is ubiquitous. Such tests can be attempted us- 
ing several approaches. For example, a particular alignment 
hypothesis can be tested with weak lensing measurements. An- 
other statistical test can be done using the projected distribu- 
tions of satellite galaxies in large galaxy surveys. In the Ap- 
pendix, we present the distribution of expected satellite angles 
with respect to the major axis of the projected host halo ellipse 
(a more complete study will be presented in a forthcoming pa- 
per, Zentner et al., in preparation). This predicted distribution 
can be compared to the observed satellite distribution around 
the major axis of the light distribution of disk galaxies. Such 
a comparison can test any correlation between disk orientation 
and dark matter elongation. A detection would indirectly in- 
dicate a preferential orientation of disk galaxies in DM halos, 
while a null result would have the interesting implication that 
the correlation between disk and halo orientation is, at most, 
weak. These tests should be feasible with existing large galaxy 
surveys, such as the Sloan Digital Sky Survey. We therefore 
expect progress in this direction in the near future. 
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APPENDIX 
SATELLITE ANISOTROPY IN PROJECTION 

Above, we demonstrated that satellite halos are not dis- 
tributed isotropically about their main host halos, but are pref- 
erentially aligned along the major axes of their hosts. This, 
coupled with the fact that the luminous dwarfs likely form in a 
biased subset of DM subhalos, greatly mitigates any claims that 
the CDM paradigm of structure formation is in conflict with the 
observed MW satellites. Agreement seems to require that the 
rotation axes of disk galaxies are nearly aligned with the major 
axes of their host halos. It would be interesting to examine pop- 
ulations of satellites around other galaxies and to compare them 
to theoretical predictions. This requires examining the angular 
distributions of satellites in two-dimensional (2D) projection, 
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where the anisotropy is somewhat diluted. We present the pro- 
jected distribution of subhalos relative to their host halos in this 
Appendix. 

To determine the projected subhalo angular distributions, we 
take three orthogonal projections of each host halo and find the 
principal axes of the 2D mass distributions by diagonalizing the 
2D analog of the moment of inertia tensor [Eq. (1)]. For each 
satellite, we determine the angle <j>, between the major axis of 
inertia and the satellite position. In practice, correlated galaxies 
along the line-of-sight to the host galaxy contribute to the differ- 
ential angular distribution as it is difficult to remove objects at a 
fixed three-dimensional distance. To account approximately for 
this, we include in our sample all halos and subhalos within a 
three-dimensional distance r cut = 2R\%q from the host halo cen- 
ter, and within a projected, 2D distance r p = 300/z~'kpc. It is not 
possible to extend our projection to larger three-dimensional 
radius because 2R^q marks the maximum extent of the high- 
resolution regions of the simulation (see § 2.1). However, in 
a forthcoming study we examine the angular distributions of 
halos relative to nearby hosts in projection and find that the dif- 
ferential distributions are not strongly affected by extending the 
projection region (Zentner et al., in preparation). 

Figure A9 shows differential distributions of <j> for the satel- 
lites of our simulated halos. The anisotropy in 2D projection 
remains clear. The probability that a satellite lies at the major 
axis (<f> < 10°) of the ellipse is > 50% larger than the probabil- 
ity that it lies near the minor axis (<f> > 80°). Additionally, the 
degree of anisotropy appears to be weakly dependent on sub- 
halo size over the relatively small range of subhalo V max that we 
can reliably probe. 

The coordinate system used in any observational study will 
be defined by the distribution of light rather than by the un- 
seen distribution of DM. As such, Fig. A9 requires a strong 
correlation between the orientation of the luminous component 
of the host galaxy and the principal axes of its host halo. In 
the absence of such a correlation, any anisotropy would be di- 
luted away in a study of many host systems. This fact provides 
an interesting constraint on theories of galaxy formation, as it 
can potentially serve to provide some statistical measure of the 
orientations of disk and/or elliptical galaxies within their host 
halos. 
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